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Rudder  angle 

deg,  rad 

P 

Air  density 

slug/ft^ 

y 

Flight  path  angle 
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Abst  rar  t 


quant  il.it  i  ve  measure  of  linear  system  controllability  quality 
was  developed,  implemented  In  a  software  program,  and  demonstrated 
through  a  preliminary  evaluation  of  an  experimental  aircraft  flight 
control  svsiem.  Aerodynamic  data  (stability  and  control  derivatives) 
st  calculated  with  linear  aerodynamics  software,  then  used  in 
the  linearized  perturbation  motion  equations  to  generate  a  linear  state 
space  representation  of  aircraft  motion.  This  representation  was  then 
used  as  the  differential  constraint  equation  for  a  quadratic  cost 
functional  which  minimizes  the  control  energy  required  to  reach  a  specified 
final  state. 

The  singular  value  decomposition  of  the  controllability  Grammian, 
which  appears  In  the  solution  to  this  optimal  control  problem,  yields  a 
unique,  time-varying,  orthogonal  set  of  basis  vectors  (Vf^,  iml,2,...,n) 
which  are  ranked  by  the  control  energy  required  to  move  the  linear  sys¬ 
tem  along  each  vector.  Eachlf^  is  the  least  costly  (in  control  energy) 
direction  to  control  in  the  orthogonal  complement  to  that  subspace  of 
(R°)  containing  \Tj ,  J»l, 2 , . . . , (i-1) . 

The  experimental  flight  control  system  consists  of  variable  inci¬ 
dence  wingtips  acting  as  elevens  and  rudders,  used  on  a  supersonic- 
cruise  lightweight  fighter.  The  evaluation  included  both  comparison  of 
basis  vector  sets  between  baseline  and  modified  aircraft  and  direct 
comparison  of  the  control  magnitudes  required  to  reach  sample  final 
states.  Results  show  the  experimental  controls  are  unsuitable  due  to 

1)  poor  longitudinal  and  lateral  control  power 

2)  no  yaw/rudder  control  available  due  to  symmetry  restrictions 
In  the  linear  aerodynamics  software. 

xv 


CALCULATING  AIRCRAFT 
CONTROLLABILITY  QUALITY 


I .  Tnt  roduct  ton 


Back^  rmi  ud 

In  the  field  of  linear  system  analysis,  the  concept  of  controlla¬ 
bility  is  well  known;  i.e.,  it  is  determined  by  the  rank  of  the  con¬ 
trollability  matrix 


;  A"" 8j 


(1) 


of  the  linear  system  described  by  the  state  space  matrix  equation 


X  -  Ax  +  Bu 


(2) 


However,  this  measure  of  linear  system  behavior  is  largely  qualitative. 
Either  the  controllability  matrix  is  of  full  rank  or  it  is  not;  there¬ 
fore,  the  linear  system  is  either  completely  controllable  or  it  is 
not.  Although  there  is  a  transformation  matrix  which  separates  the 
state  space  into  completely  controllable  and  completely  uncontrollable 
modes,  this  is  still  only  a  qualitative  result.  The  only  Information 
gained  about  the  system's  controllability  characteristics  is  a 
"yes-or-no"  answer  as  to  whether  a  given  mode  la  controllable  or  not. 

In  terms  of  "ease"  or  "difficulty"  of  controlling  the  controllable 
modes,  this  method  gives  no  quantitative  measure  of  system  controlla¬ 
bility.  Neither  is  there  any  quantitative  indication  of  "how  close  to 
uncontrollable"  are  the  controllable  modes. 

Compare  this  situation  with  that  of  stability  analysis.  The 


1 


techniques  used  in  stability  analysis  not  only  determine  the  qualita¬ 
tive  "yes-or-no"  answer  concerning  stability  of  each  mode,  but  also 
quantitatively  describe  the  characteristic  modes  of  system  motion 
(e.g.,  exponential  growth,  exponential  decay,  damped  sinusoid,  etc.). 
In  addition,  there  are  quantitative  methods  to  measure  the  stability 
characteristics  of  the  system,  such  as  time  constants,  damping  ratio, 
etc.  We  now  explore  a  similar  quantitative  measure  of  system  con¬ 
trollability  characteristics  which  would  be  useful  as  a  control 
system  analysis  tool. 

This  study  developed  a  method  to  identify  that  unique,  time- 
varying  set  of  characteristic  directions  in  state  space  (Rn)  which 
describes  the  system  controllability  quality.  In  this  set,  the  n 
directions  are  mutually  orthogonal  and  therefore  the  set  spans  (Rn) . 
Also,  these  directions  can  be  ranked  according  to  the  relative  amount 
of  control  energy  required  to  reach  a  fixed-magnitude  final  state 
located  along  each  respective  direction.  Not  only  can  theae 
directions  be  ranked  with  respect  to  each  other,  but  they  also  indi¬ 
cate  the  easiest  and  hardest  to  attain  directions  in  the  entire 
system  state  space  (Rn).  The  definitions  of  these  special  directions 
are: 

Direction  1)  the  easiest  direction  of  all  to  attain  (i.e., 
least  control  energy  to  reach)  in  the  entire  system  state  apace  (Rn) ; 

Direction  2)  the  easiest  direction  of  all  to  control  (least 
control  energy  to  reach)  in  (R°-*) ,  that  particular  subapace  of  (Rn) 
which  is  the  orthogonal  complement  to  the  one-dimensional  subspace 
spanned  by  direction  1. 
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Direction  i)  the  easiest  direction  of  all  to  control  in  (Rn 
that  particular  subspace  of  (Rn)  which  is  the  orthogonal  complement 
to  the  (i-l)-dimensional  subspace  spanned  by  directions  1  through 

(i-D. 


Direction  n)  the  most  difficult  direction  of  all  to  control 
(most  control  energy  to  reach)  in  the  entire  state  space  system  (Rn) . 
This  new  quantitative  measure  of  controllability  quality,  then,  con¬ 
sists  of  the  relative  amounts  of  control  energy  required  to  move  the 
linear  system  along  each  orthogonal  direction  in  this  unique,  time- 
varying  set. 

The  development  of  this  quantitative  measure  of  linear  system 
controllability  appears  in  detail  in  Chapter  IV  and  employs  the  follow 
ing  rationale.  First,  consider  the  functional 

J(ult))  5  *  ir)  <lr  (3) 

o 

When  the  matrix  R  is  limited  to  diagonal  form,  with  only  non-negative 

-T  -  *  2 

diagonal  elements,  the  Integrand  (u  Ru)  becomes >■  u  R. which  may 

be  interpreted  physically  as  "control  energy",  analogous  to  the 

2 

expression  for  kinetic  energy,  >»  mV  .  Restricting  R  to  non-negative 

diagonal  form,  then,  the  entire  functional  expression  J(u(t)>  is  the 

net  control  energy  used  during  the  time  interval  [0,t^].  Now  let  the 

Linear  system  begin  at  zero  initial  state  and  alternately  drive  it  to 

the  same  final  state,  x{,  at  the  same  final  time,  tf,  by  using  two 

-1  -2 

different  control  Input  histories,  u  (t)  and  u  (t),  which  result  in 
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-1  -2 

two  different  respective  state  space  trajectories,  x  (t)  and  x  (t). 

In  general,  one  of  these  control  histories  will  require  less  control 

energy  than  the  other,  as  measured  by  the  relative  magnitudes  of 

-1  -2 

their  respective  functionals,  Jj*J(u  (t))  and  l^JCu  (t));  i.e.,  for 
a  fixed  final  time  and  final  state,  the  control  history-response 
trajectory  combination  with  the  smaller  J(u(t))  is  "easier"  to  attain 
since  it  requires  less  control  energy. 

Second,  we  use  this  control  energy  functional  to  make  the  quanti¬ 
tative  evaluations  of  the  controllability  characteristics  of  the 
linear  system.  Rather  than  comparing  arbitrary  control  histories 
which  drive  the  system  to  an  arbitrary  final  state,  we  wish  to  select 
those  specific  final  states  which  can  be  used  to  identify  absolute 
limits  of  control  system  capabilities.  The  minimum  control  energy 
is  found  by  solving  the  following  constrained  optimization  problem: 
Find  the  u(t)  that  minimizes  the  quadratic  functional 

if 

s  i  J  a  fcu  ctj  <Jr  (4) 

e 

where  R  i6  positive  definite,  subject  to  the  linear  differential  con¬ 
straint 

=  A  ■*  +  B  a 

where  x(0)  =  0  and  xf  •  x(t^)  are  specified.  Using  the  singular 
value  decomposition  of  the  controllability  Graamian, 

*  r 

vtc  (t)  s  j  B8  # m  <ir  (5) 

0 

where  ^  (t)  is  the  state  transition  matrix,  e^C,  we  calculate  the 
unique  set  of  n  characteristic  directions,  defined  above,  in  (Rn) 
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which  describe  system  controllability  quality. 

Application  of  the  Concept 

To  demonstrate  one  possible  application  of  this  quantitative  con¬ 
trollability  measure,  consider  a  "controllability  comparison"  between 
two  different  aircraft.  Presently,  the  qualitative  determination  that 
one  aircraft  is  controllable  while  the  other  is  not  is  the  only  rela¬ 
tive  measure  of  controllability  performance.  If  the  two  aircraft  are 
either  both  controllable  or  both  uncontrollable,  no  information  has 
been  gained  from  which  to  decide  any  relative  merit  of  either  aircraft 
system.  In  such  cases,  we  have  no  quantitative  information  as  to 
which  aircraft  is  "more  controllable".  A  typical  application  of  this 
new  quantitative  controllability  measure,  then,  is  to  generate  these 
quantitative  results  for  different  aircraft  and  compare  them  in  order 
to  deduce  some  general  conclusions  about  relative  controllability 
performance  between  the  two  aircraft-control  system  combinations. 
Rather  than  compare  two  totally  different  aircraft,  this  study  con¬ 
sidered  the  comparison  of  two  different  control  systems  Installed  on 
the  same  basic  aircraft.  Specifically,  this  was  a  preliminary  eval¬ 
uation  of  a  new  flight  control  system  concept  compared  to  a  more 
conventional  flight  control  system. 

Control  System  Concept  and  Aircraft  Selection 

Most  current  aircraft  flight  control  systems  involve  three 
separate  sets  of  control  surfaces,  each  set  primarily  producing 
moments  about  only  a  single  axis  of  rotation  of  the  eircraft.  An 
alternative  to  this  one-axis  moment  control  surface  is  the  all-flying 
horizontal  tail  which  generates  both  rolling  and  pitching  moments  by 
using  asymmetric  and  symmetric  stabilizer  deflections,  respectively. 
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Logically,  the  simplest  moment  control  system  using  such  a  multi-axis 
concept  would  be  a  single  surface,  or  set  of  surfaces,  which  would 
provide  moments  about  all  three  axes  of  rotation.  The  experimental 
control  system  tested  in  this  study  is  an  attempt  to  use  such  a 
system,  and  Is  configured  as  follows. 

The  new  flight  control  surfaces  are  the  wlngtlps  (Figure  1), 
which  act  as  conventional  elevons  to  generate  pitch  and  roll  control, 
and  are  deflected  individually  to  produce  yawing  moment.  Two  different 
forces  produce  yawing  moments  in  this  case:  first,  asymmetric  drag  of 
the  single  deflected  tip;  and  second,  side  force  of  the  deflected  tip. 

Using  the  moment  equation 

M  -  R  x  F  +  R  x  F  (6) 

z  y  x  x  y 

Figure  1  shows  the  wLngtip  drag  force,  acting  through  a  y-axis 


Figure  1.  Yawing  Moment  Due  to  Individual  Tip  Deflection 
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moment  uni,  together  wltli  Fj_  ,  the  force  component  of  the  net 

wingtip  fill  force  l'L  ,  acting  through  an  x-axls  moment  arm,  combining 
to  prodint*  the  resultant  M  yawing  moment. 

The  aircraft  chosen  fur  this  new  control  system  is  the  Boeing 
Lightweight  Experimental  Supercruiser  (LES)  model  230-Y  (Figures  2 
and  3).  There  are  two  reasons  for  this  choice.  First,  geometry  data 

was  available  from  Boeing  for  modeling  the  230-Y  with  FLEXSTAB  to  calc¬ 
ulate  ahl  necessary  aerodynamic  data.  Second,  as  a  contractor-funded 
design,  the  320-Y  represents  a  practical,  physically  realistic  baseline 
control  system  against  which  to  compare  the  new  system. 

Purpose 

The  purpose  of  this  study  was  to  both  develop  a  quantitative 
measure  of  linear  system  controllability  characteristics  and  to 
demonstrate  this  measure  as  the  main  tool  in  a  preliminary  evaluation 
of  an  unorthodox  aircraft  flight  control  system. 

Coals 

There  are  three  separate  goals  of  this  study.  The  first  goal  is 
the  theoretical  development  of  this  quantitative  measure  of  linear 
system  controllability.  The  second  goal  is  the  implementation  of  the 
software  to  calculate  this  quantitative  measure.  The  third  goal  is 
the  practical  demonstration  of  this  scftvara  by  both  using  it  to  cal¬ 
culate  the  degree  of  controllability  of  three  different  aircraft/ 
flight  condition  combinations,  then  usi'.ng  these  results  to  make  a 
preliminary  evaluation  of  the  particular  flight  control  systems  used 
on  the  experimental  aircraft. 

Assumptions 

First,  this  study  assumed  a  rigid  aircraft  ignoring  the  effect  of 
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structural  dynamics  on  the  design  of  the  control  systems.  The  second 
restriction  was  to  address  only  linear  systems,  which  permits 
linearization  of  both  dynamics  and  aerodynamics.  The  fairly  uncon¬ 
ventional  aircraft  configuration  chosen  for  analysis  reasonably 
precluded  the  use  of  the  theoretical/empirical  DATCOM  methods 
developed  from  more  orthodox  configurations.  Instead,  the  FLEXSTAB 
(References  1  and  2)  digital  computer  program,  based  on  linear  aero¬ 
dynamic  theory,  was  used  in  this  study.  Although  wind  tunnel  testing 
would  provide  the  most  accurate  aerodynamic  data,  the  linear  modeling 
of  aerodynamics  should  not  significantly  impair  the  accuracy  of  the 
results,  since  the  study  Is  limited  to  the  same  small  perturbation 
region  about  some  equilibrium. 

Approach  and  Overview 

The  first  step  in  the  study  was  to  use  the  FLEXSTAB  program  to 
estimate  the  aircraft  stability  and  control  characteristics.  This  data 
is  presented  and  discussed  in  Chapter  II. 

The  second  step  was  to  convert  the  stability  and  control  deriva¬ 
tives  from  FLEXSTAB,  using  the  linearized  equations  of  motion,  into 
the  constant-coefficient  matrices  of  the  linear  system  state  space 
representation.  The  results,  along  with  a  derivation  of  the  linearized 
aircraft  equations  of  motion,  are  given  in  Chapter  III. 

The  third  step  was  to  quantitatively  measure  the  degree  of  con¬ 
trol  capability  of  each  linear  system  at  the  defined  flight  condition. 
Theoretical  development  of  this  quantitative  measure  Is  presented  in 
Chapter  IV.  The  results  of  the  control  system  evaluation,  in  terms 
of  this  measure  of  controllability,  are  given  In  Chapter  V.  Chapter  VI 
lists  the  conclusions  and  recommendations  of  this  study. 
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1  l .  Aerodynaml  r  Da_t. a  Genera t  Ion 
Basic  Theory  and  Outputs 

Aircraft  flight  character isties  were  calculated  using  the  FLEXSTAB 
computer  program  (References  1  and  2).  FLEXSTAB  uses  an  aerodynamic 
theory  based  on  a  first-order,  linear,  small-perturbation  approximation 
to  the  unsteady,  inviscid,  irrotut ional,  velocity  potential  equation: 

where  a  is  the  local  speed  of  sound  for  both  subsonic  and  supersonic 
flow.  This  partial  differential  equation  has  the  physical  boundary 
condition  of  zero  normal  velocity 

V-VF  -  0  (8) 

at  F(x,y,*,t)  •  0  ,  the  aircraft  surface. 

Details  of  this  method  are  presented  in  Reference  1. 

The  net  FLEXSTAB  outputs  are:  (1)  for  a  straight  and  level  flight 
condition  at  the  trim  airspeed  U^,  with  all  other  trim  parameters  and 
control  surface  deflections  set  to  zero,  the  aerodynamic  force 
derivatives  are  calculated;  (2)  these  derivatives  are  used  in  linear¬ 
ized  equilibrium  equations  to  calculate  the  actual  trim  values  for  all 
the  aircraft  state  and  control  variables;  and  (3)  these  trim  values  are 
finally  used  to  calculate  the  stability  and  control  derivatives  for  the 
trimmed  reference  flight  condition. 

Control  Surface  Definitions 

FLEXSTAB  assumes  the  aircraft  to  be  symmetrical  with  respect  to 
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the  XX  plane,  as  in  Figure  4,  and  only  accepts  control  surfaces 


located  either  on  the  XX  plane  or  else  located  symmetrically  with 
respect  to  the  XZ  plane.  Also,  control  surface  pairs  not  on  the  XZ 
plane  must  be  deflected  symmetrically,  like  elevators,  or  asymmetri¬ 
cally  ,  like  ailerons.  The  effect  of  a  single  surface  deflection  on 
either  side  of  the  XZ  plane  cannot  be  modeled.  This  PLEXSTAB  restric¬ 
tion,  then,  reduces  the  experimental  flight  control  system  to  a  pair 
of  wingtips  acting  as  elevons,  with  no  provision  for  rudder  effects 
except  for  the  coincidental  lateral/directional  effects  of  ordinary 
aileron  deflections.  For  conventional  control  systems,  the  following 
sign  convention  for  control  deflections  is  adopted:  down  elevator, 
right  aileron  up,  and  left  rudder  deflections  are  taken  as  positive. 
For  the  experimental  control  system,  vingtip  deflections  and 

£  t are  defined  as  positive  for  vingtip  trailing  edge  down. 

Applying  this  convention  to  all  three  moment  axes,  vingtip  deflections 
generate  the  same  moments  as  the  following  conventional  controls: 

positive  &  tipR  ectu  ivaleoC  to  positive  negative  £A 

positive  £  equivalent  to  positive  positive 


Using  the  following  definitions  for  control  surface  deflections 
for  the  unconventional  aircraft. 


t( 

&  r  + 

£  T>*  ) 

4. 

(9) 

ft 

< 

^0 

£r,<v  ) 

(10) 

M 

+ 

)  •  'I* 

(11) 

it  is  apparent  that  there  is  no  "pure"  rudder  deflections  possible. 
Any  roll  or  yaw  effects  due  to  rudder  are  already  contained  in  the 
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TIP  (THIN  BODY  #2) 


WING  (THIN  BODY  #1) 


FUSELAGE 
(SLENDER  BODY) 


Figure  4.  Aerodynamic  Paneling  Approximation 


aileron  derivatives,  and  any  pitch  el  lefts  due  to  rudder  are  already 
eontained  in  the  elevator  derivatives. 

Since  some  rudder  surface  must  be  defined  in  FLEXSTAB,  a  strak*  was 
included  as  an  artificial  rudder  to  allow  FLEXSTAB  to  always  successfully 
trim  the  experimentally-controlled  aircraft,  but  all  strake-related 
rudder  control  derivatives  were  discarded. 

Discussion  of  Data 


Tables  A.T,  A. II,  and  A. Ill  show  the  respective  FLEXSTAB-computsd 
results  tor  both  the  baseline  and  experimental  LES  230-Y  aircraft  at 
the  indicated  flight  condition.  Comparing  the  two  LES  230-Y  models, 
the  differences  in  several  trim  and  stability  parameters  warrant 
explanation:  (1)  trim  values  for  angle  of  attack,  thrust,  and  eleva¬ 

tor;  and  (2)  static  longitudinal  stability,  C  .  The  trim  angle  of 

“•l 

attack  difference  is  probably  due  to  the  baseline  aircraft's  wingtip- 
mounted  vertical  fins.  They  provide  an  endplate  effect  which  generates 
a  stronger  lift  distribution  than  for  the  fin-less  experimental 
aircraft.  The  baseline  aircraft  can  thus  produce  the  same  amount  of 
lift  at  a  smaller  angle  of  attack.  The  baseline  aircraft's  larger 
trim  thrust  is  also  due  to  its  vertical  fins,  which  produce  more  drag 
(thereby  requiring  more  thrust)  than  the  fin-less  experimental 
aircraft.  This  erroneously  large  disparity  in  trim  thrusts  is  due  to 
the  inherent  inability  of  FLEXSTAB  linear  theory  to  accurately  predict 
drag  data.  Experimental  wingtlp  trim  deflection  of  35s  would  certainly 
produce  even  more  drag  and  require  more  trim  thrust  than  would  the 
baseline  fins  plus  elevators  trimmed  at  only  48  deflection. 

The  difference  in  trim  elevator  settings  is  due  to  the  differences 


between  the  two  types  of  elevator  control  surfaces.  Table  A. Ill  shows 
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iu.it  itu  i-xpii  inu-ut  .1 1  I'lfvuiit;  jic  I  .ir  loss  effective,  as  measured  by 
(.'  and  C,  .  flie  reasons  lor  Lliis  include:  (a)  the  experimental 

“r  Lr 

&6  "E 

wingtip  elevons  are  physically  much  smaller  In  surface  area  than  the 
baseline  trailing  edge  elevons,  so  can  generate  less  change  in  net 
lift  pot  degree  of  deflection;  (b)  the  control  surface  center  of 
pressure  for  the  baseline  elevons  has  a  longer  moment  arm,  with  respect 
to  the  aircraft  center  of  gravity,  than  for  the  experimental  wingtip 
elevons,  so  the  experimental  surfaces  produce  less  change  in  pitching 
moment  for  the  same  amount  of  change  in  net  lift;  (c)  the  extreme 
wingtip  anhedral  (60°)  reduces  by  half  the  vertical  component  of  the 
change  in  net  lift  used  to  generate  changes  in  pitching  moment.  These 
three  effects  combine  to  severely  limit  the  experimental  wingtip 
elevons'  pitch  effectiveness.  Since  an  aircraft  would  never  be 
designed  to  fly  with  such  large  trim  settings,  a  better  preliminary 
design  approach  would  be  to  modify  the  trailing  edge  camber  of  the 
wing  (i.e.,  to  "design  in"  more  up-elevator)  to  allow  a  more  reason¬ 
able  trim  elevator  setting  for  the  experimental  aircraft.  Also,  more 
control  power  could  be  gained  by  increasing  the  physical  size  of  the 
wingtip  elevons. 

The  experimental  aircraft  is  also  statically  unstable  in  pitch, 
as  seen  by  a  positive  C  in  Table  A. III.  Both  LSS  models  vara  balancad 
at  the  same  percent  mean  aerodynamic  chord,  and  in  fact  have  the  same 
identical  wing  geometry,  M.A.C.,  and  longitudinal  location  of  M.A.C. 
FLEXSTAB  output  data  shows  chat  the  trim  neutral  point  was  slightly 
behind  the  CG  for  the  baselLne  aircraft  (stable),  and  was  slightly 
ahead  of  the  CG  (unstable)  for  the  experimental  aircraft.  This 
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neutral  point  sliilt  between  identically-winged  aircraft  is  also  due  to 
tlve  presence  of  vertical  fins.  Since  tULs  study  was  not  directed 
toward  aircraft  design,  these  results  were  left  as  is,  and  the  analysis 
continued  after  noting  the  problems  discussed  above. 

The  stability  and  control  derivative  data  of  table  it. Ill  is 
now  converted  to  the  linear  state  space  representation  of  aircraft 
perturbat ion  motion  by  the  conversion  process  described  in  the  next 
chapter. 
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lonvcrtlng  Fl.F.XSTAB- gene  rated  stability  and  control  derivatives 
Into  the  linear  state  space  representation  for  the  various  aircraft- 
flight  condition  combinations  begins  with  the  six  nonlinear  differen¬ 
tial  equations  of  motion  and  Che  two  kinematic  relationships : 
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The  first  step  Is  to  form  the  perturbed  state  equations  of 
motion  with  respect  to  steady  state  flight.  Steady  state  flight 
Implies  (14:2.25)  that  all  motion  variables  remain  constant.  The 
perturbation  substitution  (for  each  force  and  moment  component  and  for 
all  the  motion  variables)  expresses  the  instantaneous  value  of  each 
variable  as  the  sum  of  the  equilibrium  value  (i.e.,  trim  value) 
denoted  by  the  subscript  (  ) ^,  plus  the  perturbation  away  from  the 

equilibrium  value  (denoted  by  the  small  case  symbol).  Table  1  shows 
the  complete  list  of  perturbation  substitutions  for  all  parameters 
Involved. 

The  second  step  is  to  apply  the  two  fundamental  linearity  re¬ 
strictions  on  the  magnitude  of  all  perturbations:  (1)  the  small  angle 
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(2)  all  perturbation  magnitudes  are  small  enough  that  any  produce  or 
cross-product  of  perturbations  is  negligible  with  respect  to  the 
perturbations  themselves.  Removing  the  equilibrium  terms  and  neglect¬ 
ing  all  nonlinear  perturbation  terms  produces  the  following  linearized 
equations  of  perturbed  motion  in  stability  axes: 
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♦  C  )  A  *  (C  “  )  +  (c«  k  )  -p  +  (X  )  V-  +  C.  *  ♦  ^  1  (26) 

+  l-Hr  JP  ^  ^ut  1  V  v  *f  aut  '  T  v  «„  ou#'  4^  -1 


(23) 


(24) 


(25) 
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s 


(27) 


^  Coi  $  -  Y  I,  -  (Uji  5C^  ^  ^  C<>S  $t)  + 

-f  +  (%  <"*  $t*  rto*  /j  +  ($1  C*'  ft~  #t  )^J  + 

t  ^l^i’i)  &  34C1  *at  (28) 

£  ^  (-Ut  ♦(U1uv^1)s  (29) 


The  third  step  is  to  substitute  the  FI.EXSTAB-generated  stability 
and  control  derivatives  into  these  linearized  equations  of  perturbed 
motion  (21)  -  (29)  as  the  constant  coefficients  of  the  motion  and 
control  variables 


T 

((*»*>  f>r)T 

(30) 

(31) 

Solving  this  system  for  the  time  derivatives  of  the  state  variables 
x  yields  the  state  space  representation  of  the  aircraft  perturbation 
motion  with  respect  to  steady  state  conditions: 


x  ■  Ax  +  Bu 


(32) 


where  the  coefficient  matrices  A  and  B  are  made  up  of  combinations  of 
the  constant  coefficients  from  the  original  motion  and  kinematic 
equations. 

Tables  11  and  III  show  the  resulting  A  and  B  matrices  calculated  for 
both  LES  230-Y  models  tested,  for  both  the  longitudinal  and  lateral/ 
directional  modes.  These  A  and  B  matrix  pairs  define  the  physical 
differential  constraint  equation  of  the  optimal  control  problem  dis¬ 
cussed  in  the  next  chapter.  All  longitudinal  and  lateral/directional 
mode  eigenvalues  of  both  aircraft  A  matrices  are  included  in  Table  IV. 
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Table  II  Baseline  Boeing  LES  230-Y  State  Space  Matrices 


\ONG  *  (£e<DEG)  »eT(*  a^craft  weight)) 

“lat  *  f  <1^( DEG)  ,^(deg)]t 


Table  IV 

Eigenvalues  of  Aircraft  A  Matrices 


AIRCRAFT 

longitudinal 

EIGENVALUES 

LATERAL-DIRECT 10NAL 
EIGENVALUES 

0.0 

Baseline 

-0.016 

.005526 

LEE 

0.0289 

-0.1602  +  j  2.307 

-0.4326  +  j  1.2 

-0.8415 

Experimental 

1.01 

LES 

-.02829  +  j  .1118 

1.119 

0.0 

-.88  +  j  0.5634 

-1.733 

.0006578 
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IV.  1'onlrol  S\  ■ .  t  fin  Analysis 


l  lii'  Might  .  out  rol  systoms  art-  now  evaluated  quantitatively  by 
their  ability  to  maneuver  the  aircrutt  from  the  reference  flight  condi¬ 
tion  to  a  spec i tied  1 inal  flight  condition  differing  from  the  refer¬ 
ence  by  nitiail  perturbations  in  the  motion  variables  x  of  equation  vj2). 
lot  a  given  reference  I  light  condition,  equation  (32)  models  the 
physical  changes  in  aircraft  motion  x  caused  by  control  deflections  u, 
when  the  instantaneous  flight  condition  is  sufficiently  close  to  the 
reference  condition  characterized  by  equation  (32).  The  quantitative 
measure  developed  in  this  study  Is  discussed  in  the  optimal  control 
technique  described  below 
Ba sic  Theory 

The  problem  is  to  select  the  optimal  control  time  history,  u(t), 
which  minimizes  the  functional  performance  index 

*  V  T 

Jiiiti):  i  1  u(t)RU(*)Jtr  03) 

where  R  is  a  positive  definite  matrix,  and  where  t^  is  specified  and 
both  x(t)  and  u(t)  are  subject  to  the  constraint  equation, 

x  «  Ax  +  Bu,  (34) 

and  x(tQ)  *  Xq  =  0  and  x(tf)  -  are  specified.  Since  x(t)  can  be 
expressed  as 

x(t)  =  $(t)x  +  f  $(t-T)Bu(?)d?,  (35) 

0  o 

the  state  variables'  time  history,  x(t),  is  dependent  on  the  control 
variablts'  time  history,  u(t).  We  can  treat  them  as  independent,  for 
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tin-  (iiii'I'om'  ul  ini;  i  In.-  I  irsL  iii-.i-ssjiy  condition  tor  optimality, 

by  introducing  n  Lagrange  multipliers^  to  form  the  augmented  perfor- 


iiaiu'e  index 


re  K 


TA(u(t))~  i  JiUr)  £urr|  it  +  ify-V  Birr)  Jr -jf  ?  (36) 

o  a  f 

Since  the  const  taint  is  introduced  involving  Lhe  constant  vector 
x  ,  >  is  also  a  constant  vector.  Combining  the  separate  integrands  of 
equation  (76)  yieLds 

T 

-•  i  J  (  u  It,  Ru  ft,  +  X  $  (tf-T)  -K  %  (37) 

Defining  the  Hamiltonian  function  as 

—  |  _  . T  __r  r 

%  (tin,  A)  *2  U.<t)Ru(t)  +  * Bu  171  (38) 


we  now  express  the  augmented  performance  index  as 


r  df  -  'K  ■% 

e  • 


The  superscript  (  )*  is  hereafter  used  to  denote  quantities  evaluated 
at  the  optimal  condition,  i.e.,  at  the  minimum  value  of  JA(u(t)). 

Examining  a  non-optimal,  neighboring  control  u(t),  lying  close  to 
the  optimal  u*(t),  but  differing  from  the  optimum  by  a  variation 
£  u(t) , 

u(t)  -  U  (11  f  Id(t)  =  U  rt)  ♦  e  (40) 

where  l(j  is  a  small  number.  Note  that  (t)  may  be  chosen  arbitrarily. 
Note  also  that  both  the  initial  and  final  conditions  x^  »  0  and 
x(tj.)  =  Xj.  are  always  optimal,  since  they  are  located  on  the  optimum 
trajectory  x*(t)  connecting  x^  and  x^..  From  the  calculus  of 
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v.u  i.n  ions,  jiiv  v.iliu-  ot  J  (u)  ,  where  n  t  u* ,  must  he  greater  than  or 

equal  to  the  optimum  (minimum)  J  (u)*.  Tite  necessary  conditions  for 

optimality  ot  u*(t)  are  found  by  expanding  J  (u)  in  a  Taylor  series 

A 

expansion  about  the  optimum  J  ( u* ) , 

1 ; 

JA«0=  J  y.[(u\^),  *]  Ar  -  ^Tx  (ad 

o 

Lxpjiul  i  ng  tint  her,  we  have 

/  *P  DtJJ  .  T 

Ja(51^  U*t-CS',Xi>  5f/  ef(t))  + 

©  »  ' 

“  (42) 
The  underlined  terms  are  just  JA(u*) ,  and  the  €  is  factored  out  to  give 

-Va*)  f[  +  b^1)  (43) 

2 

Since  €.  can  be  chosen  arbitrarily  small,  all  terms  ot  order  €  and 
higher  can  be  made  negligibly  small  with  respect  to  the  €  -term.  Also, 
since  <  can  be  chosen  Co  be  either  positive  or  negative,  the  only  way 
JA(u)  can  always  be  greater  than  or  equal  to  JA(u*) ,  for  all  permissible 
(  and  ">[(t),  is  to  require  that  the  coefficient  of  the  #  -term  in 
equation  (43)  always  be  zero.  This  necessary  condition  for  u*(t)  to 
be  the  optimal  control  is  then 


J  J*  ' 
o  * 


Applying  the  fundamental  Lemma  of  the  calculus  of  variations,  the 
necessary  condition  for  optimality  becomes 


x)j 

u*(»  =  -R-1BT^T(tf->)  > 


=  0,  or 
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I'  .i  mm!  in-'  r>;,,.ii  i, hi  (4t  I  ini  h(i  >  Inin  eouation  (3j)  gives 


*«  T 

*(tf)  =  -  o  $<**'*)  Br’ Br  $(i<-r) 

T  o 


(47) 


l»\  i  lui  i. .  i  ii|',  i  In  dummy  vjt  iable  of  Liu  igrai  ion  lroin  to  (t^-t),  the 
mt  m •  •  *  I  in  i  iiu.it  ion  (47)  becomes 


J  $(n  8R*8 


Kspri'st,  ing  the  comb  i  nation  BH  *Bl  as  B  B  where 

K  K 

R«  =  br-j 


(48) 


the  integral  becomes 

T  r 

r  ^  "c 

o 

This  integral-matrix  expression  is  defined  in  the  literature  (8:402) 
as  the  symmetric  (n  x  n)  controllability  Crammian,  W^(tj ) .  For  both 
t^  =  0  and  t  *  tj.  specified  in  this  optimization  problem,  the  bracketed 
integral  in  equation  (47)  becomes  W^(tj.).  Substituting  this  form  for 
the  integral,  equation  (47)  becomes 

xuf)  =  "VV  ^  (49) 

and  equation  (46)  becomes 

u‘(t)  =  rV§T(tf-Tr)wc“'(tf)x(tf)  (50) 

Substituting  tills  expression  for  u(t)  into  equation  (33)  gives 

t*  ,  T 

j(u(t.))  =  {^xr(tf)Wc(M  $(lf-T)8R  1  BTfT(t4-^)Wc‘1(tf)9((tf)«iT  (51) 

o 


°7 


\L  1  1  J.  lui  a 


is  involving  x(t(.)  and  W ^(tp,  art  invariant  over  the  inter¬ 
nal  oi  integration,  and  may  be  factored  out  of  the  integrand  to  give 


Thu  inu-gral  in  this  equation  is  W  (tr).  so 

C!  i 

T 

T  -i 

J(u(t))=£*4  (l*\  *  4 


(53) 


Since  W  is  a  symmetric  matrix,  the  optimal  cost  functional  is 


X<u*<0»  5 


1 

l 


*1 


(54) 


The  controllability  Grammian  can  be  decomposed  into 


wc(tf) 


pw  pw  T 


(55) 


where 


• 


is  a  diagonal  matrix  of  eigenvalues  of  W^t^),  and  P, 


H-f5! 


is  au  orthogonal  matrix  of  column  vectors  which  are  the  eigenvectors 
of  W^(t j) .  because  VJ^t^)  is  symmetric,  the  eigenvalues  and  eigen¬ 
vectors  are  the  same  as  the  singular  values  and  singular  vectors, 
respectively.  Substituting  (55)  into  (54)  gives 
♦  t  *1  T 

J  (£Ufl  =  i  Pu.  Pu,  ' 

s  i  (  VP-  )T  <56) 

-  T 

The  quantity  (x}  P^)  is  expanded  to  give 
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[< 


4  i¥. 


)  '■'  ( 


T  _ 


>J 


(57) 


which  is  a  row  vector  of  scalars,  [c  ,  c.  ,  ...  ,  c  ),  where  c.  is  the 

12  n  1 

inner  (scalar)  product  between  x^.  and  ,  i.e.  is  the  projection 
ot  x(  along  tile  eigenvector  V\  of  Wc(t^). 


CL=  <  ) 


(58) 


Substituting  this  result  into  equation  (56)  gives 


jtauiw  i[V-c-]2 


■  UK  $ 


(59) 


2 

This  is  a  weighted  sum  of  squares  (c^  )  of  the  components  of  x^.  along 

each  of  the  eigenvectors  The  individual  weights  (1/fc  ),  as 

reciprocals  of  the  eigenvalues  ^  of  W  (tr),  are  now  ordered  in 

1  c  t 

ascending  order;  i.e.,  where  ^  is  the  largest  eigenvalue  from  the 
singular  value  decomposition  of  Mc(t^),  bo  (1/A^)  is  now  the  smallest 
weighting. 

The  final  conclusion  from  this  form  of  the  cost  functional  is 
that  an  xf,  of  fixed  magnitude,  selected  along  each  of  the  "ranked" 
eigenvector  directions^!,  will  give  (in  ascending  order)  cost  values 
Jt(u*),  the  control  energy  required  to  reach  xf.  That  is,  if  xf  is  a 
unit  vector  along  V^,  then  the  control  cost  to  reach  will  be  just 
l/*l  .  If  Xj  is  a  unit  vector  along  then  the  control  cost  to  reach 
x(.  will  be  just  1/^  )  1/*^.  Also,  if  x^  is  along  <n»  the  control 
cost  to  reach  xf  will  be  l/*n  >  1/A^i  “  1,2,  ...  ,  (n-1).  This 
characterist  Ic  is  the  desired  measure  of  relative  controllability 
quality  throughout  the  entire  state  space.  That  is,  since  all 
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~f‘  ^  x  \?  j »  the  eigenvector  set  wilt  completely  spaa  the  state  space 

(Rn) ,  and  therefore  forms  a  set  of  basis  vectors  for  the  stale  space. 

Also,  sloce  the  controllability  Gramiuian  la  time-varying,  the  ranked 

eigenvector-basis  vector  Bet  la  also  time-varying.  This  implies  the 

importance  of  obtaining  time  histories  of  the  vT^(t)  vectors  over 
some  finite  time  Interval  of  interest. 

Sof tware  lmplementat Ion 

This  controllability  quality  measure,  derived  from  the  controlla¬ 
bility  Grammian,  has  been  successfully  Implemented  in  the  computer 
program  GRAMOP,  whose  output  Includes  two  sets  of  plots  (and  print): 

1)  a  set  of  n  "traits”  plots,  plus  a  cost  history  plot; 

2)  a  set  of  state  trajectory  and  optimal  control  history  plots 
for  each  user-defined  final  time  and  final  state  x(tj). 

(  The  1th  trait  plot  (e.g.,  figure  Bl)  shows  the  history,  versus  final 

a  U  aU 

time,  of  the  i  -ranked  (i  -smallest  amount  of  control  energy  re- 
l 

it  quired)  response  direction  ^(t^)  in  state  space.  Recall  from 

*  the  above  that  V^Ctj)  is  one  of  the  n  mutually  orthogonal,  time- 

varying  eigenvector/basis  vector  directions  spanning  (Rn) ,  and  that 
>T^(tj)  is  just  the  i1*1  eigenvector  of  the  controllability  Gramnian 
Wc(tj).  The  cost  history  plot  shows  the  relative  control  costs 
needed  to  reach  each  one  of  these  n  ranked  directions  in  state  space. 
The  control  and  state  trajectory  history  plots  show,  respectively: 

(1)  the  optimal  control  history  u*(t)  that  drives  the  system  from 
zero  initial  state  to  the  user-defined  final  state,  using  minimum 
control  energy  and  (2)  the  corresponding  state  trajectory  history 
x(t)  along  which  the  system  travels.  Chapter  V  presents  the  GRAMOP 
results  for  the  two  aircraft  studied. 
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V. 


l-'oi  tliis  sUuly,  two  separate  uircraf  t/f  1  ight  condition  combina¬ 
tions  were  analyzed: 

L)  LES  2  10-Y  -  Mach  0.8,  35,000  ft,  steady  level  cruise 

2)  Modi! led  1.ES  2 10-V  -  Mach  0.8,  35,000  ft,  steady  level  cruise 
The  two  l.KS  aircraft  were  modeled  using  the  FLEXSTAB  aerodynamic 
paneling  program  to  generate  stability  and  control  derivatives  for 
the  indicated  flight  conditions.  These  derivatives  were  used  to  pro¬ 
duce  a  state  space  system  of  equations  for  each  aircraft/flight 
condition  combination.  These  systems  of  equations  were  used  as  the 
constraint  equation  (34)  in  the  GRAMOP  program  in  order  to  calculate 
and  plot  the  basic  longitudinal  and  lateral/directional  controlla¬ 
bility  characteristics  of  each  aircraft  when  flying  at  the  indicated 
flight  conditions.  These  state  and  control  history  plots  are  best 
interpreted  by  first  noting  the  vertical  axis  scaling  of  each 
separate  plot  variable.  This  scaling  clearly  indicates  which  variables 
are  dominant,  in  terms  of  relative  magnitude,  in  each  response  vector 
time  history.  Each  test  case  is  discussed  in  detail. 

Aircraft  Traits 


This  Boeing-designed  aircraft  has  wing  trailing  edge  elevons  and 
rudders  on  the  wlngtip-mounted  vertical  fins,  as  shown  in  Figure  2  . 
This  aircraft  is  used  as  the  baseline  for  the  modified  LES,  to  explore 
changes  within  a  single  aircraft  design, 
a)  Direction  1  (Figure  BL) 

For  the  baseline  LES,  this  easiest -to-control  response  is  a  combination 
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i>t  piLili  i.ite,  pitch  angle,  and  angle  ol  attack  up  until  about  2 
.seconds,  when  it  becomes  almost  entirely  an  altitude  change,  with  some 
minor  velocity  change,  for  times  beyond  2  seconds.  Sample  state  and 
control  trajectories  at  1.5  seconds  (Figures  B2  and  B3)  show  very 
small  throttle  and  elevator  requirements  to  drive  the  aircraft  along 
this  direction.  1'his  is  consistent  with  a  fighter  aircraft  designed 
for  high  maneuverability. 

b)  Direction  2  (Figure  BA) 

This  response  begins  as  mostly  negative  velocity,  angle  of  attack,  and 
pitch  angle,  but  quickly  (1.0  seconds)  inverts  to  a  combination  of 
velocity  change  with  opposite  pitch  angle  and  opposite  angle  of 
attack. 

c)  Direction  3  (Figure  B5) 

This  response  exhibits  several  changes  in  behavior,  beginning  as 
mostly  pitch  angle  and  angle  of  attack,  then  (around  1.3  seconds) 
changes  to  predominantly  altitude  until  about  3  seconds,  when  it 
becomes  mostly  pitch  rate. 

d)  Direction  4  (Figure  B6) 

This  vector  begins  as  mostly  altitude,  quickly  (at  1  second)  changes 
to  mostly  velocity,  with  a  minor  and  fairly  equal-magnitude  combination 
of  pitch  angle,  pitch  rate,  and  angle  of  attack. 

e)  Direction  5  (Figure  B7) 

The  most  difficult  response  vector  remains  mostly  a  combination  of 
angle  of  attack  and  opposite  pitch  angle  through  the  final  time 
interval  tested.  The  altitude  component  dies  out  very  rapidly  (at 
about  0.5  sec.)  while  the  opposite  pitch  angle  component  grows  rapidly 
at  the  same  time.  The  difficulty  of  this  response  is  apparent  in  the 


i-x  L  riiiu- 1  \  1 .1 1  gc  eli'v.uor  control  mugn  i  t  udi's  required  to  achieve-  it 

over  shot  t  t(,  as  seen  i it  figure  B9  .  The  linearity  assumptions  of  the 
constraint  equation  are  grossly  exceeded, 
t)  Cost  Plot  (Figure  B 1 2 ) 

this  t  igure  shows  the  relative  costs  of  response  along  each  tT^(t-). 
CRAMOP  correctly  shows  that  no  two  cost  traces  overlap,  since  J(u^) 
must  always  be  less  than  j(u^)  for  j>l  at  any  final  time,  t^.  It  is 
important  to  notice  that  all  cost  lines  tend  to  decrease  sharply  as  t^ 
increases  out  to  about  2  seconds,  then  all  lines  remain  somewhat 
constant  over  later  final  times.  Also,  the  control  histories  generally 
show  progressively  smaller  control  deflections  used  as  final  time 
increases  through  the  2-second  area,  which  accounts  for  the  sharp  drop 
in  cost  magnitude  through  that  interval. 

2)  LES  230- Y  -  Lateral/Directional  Case 

a)  Direction  1  (Figure  B13) 

This  easiest  response  begins  as  almost  entirely  roll  rate,  which 
quickly  diminishes  and  is  replaced  by  roll  angle.  Control  history 
Figure  B15  shows  minimal  control  requirements,  again  attributable  to 
the  aircraft's  combat  mission  design. 

b)  Direction  2  (Figure  B16) 

This  response  direction  history  begins  as  mostly  yaw  rate,  which  dies 
out  very  rapidly  (at  about  0.3  sec.)  and  is  replaced  by  roll  angle, 
which  is  then  (at  1.5  sec.)  gradually  replaced  by  roll  rate  as  the 
dominant  component. 

c)  Direction  3  (Figure  B17) 

This  response  begins  as  mostly  roll  angle,  then  very  quickly  (at  about 
.5  sec.)  becomes  mostly  yaw  rate  and  remains  that  way  throughout  the 
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rest  ot  tin.-  t  iui«.'  interval  tested. 

d)  Direct  ion  A  (Figure  B18) 

This  response  is  mostly  sideslip  throughout  the  entire  time  interval 
tested.  Control  magnitudes  are  small,  but  usually  Increasing  steadily 
with  increasing  t  (see  Figures  B20  and  B22) . 

e)  Cost  Plot  (Figure  B23) 

The  cost  plot  characteristics  are  similar  to  the  longitudinal  case 
cost  plot . 

1 )  Modified  l.ES  230-Y  Longitudinal  Case 
The  modifications  to  the  Boeing  design  LES  230-Y  include  removal  of 
the  wingtip/vertical  fins,  elevons  "locked"  in  the  wing  surface,  and 
the  anhedral  wingtips  used  as  the  only  set  of  control  surfaces 
(Figure  3). 

a)  Direction  1  (Figure  Cl) 

This  response  begins  as  almost  entirely  pitch  rate,  then,  around  1 
second,  gradually  becomes  mostly  a  combination  of  both  some  lesser 
velocity  and  angle  of  attack  changes  along  with  a  larger  altitude 
change.  Consistent  with  this  change,  the  trend  between  the  two 
sample  trajectories,  Figures  C2  through  C5,  as  final  time  is  increased, 
is  toward  virtually  no  throttle  changes  and  only  elevator  movement. 

Note  also  that  all  control  input  histories  show  both  left  and  right 
wingtips'  movements  are  identical,  which  is  consistent  with  both  tips 
acting  together  as  elevators.  Control  magnitude  drops  with  larger 
final  times. 

b)  Direction  2  (Figure  C6) 

This  response  history  begins  as  mostly  velocity  change,  then  changes 
to  a  fairly  equal  combination  of  pitch  rate  and  opposite  altitude 
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!t  then  gradually  becomes  a  new  combi- 


change  Ik-Iuci'H  I  and  seconds, 
nation  oi  all  variables  -  positive  velocity  and  altitude  changes  along 
with  negative  angle  of  attack,  pitch  angle,  and  pitch  rate. 

c)  Direction  I  (Figure  C7) 

This  response  begins  as  a  combination  of  angle  of  attack  and  pitch 
angle,  then  changes  at  1  second  to  a  combination  of  altitude  and 
opposite  pitch  rate.  At  2  seconds  it  again  changes  to  mostly  velocity 
change  with  a  minor  pitch  rate  component. 

d)  Direction  4  (Figure  C8) 

This  response  begins  as  altitude  change,  then  (around  1  second) 
changes  to  a  fairly  constant  combination  of  pitch  rate  and  opposite 
pitch  angle,  with  minor  opposite  angle  of  attack. 

e)  Direction  5  (Figure  C9) 

This  most  difficult  response  history  is  always  mostly  angle  of  attack 
with  opposite  pitch  angle.  Figures  CIO  through  C13  show  that  this 
response  is  unattainable  within  the  linear  region  of  either  state  or 
control  motion  until  around  8  seconds  final  time. 

f)  Cost  Plot  (Figure  C14) 

Cost  plot  trends  generally  follow  those  of  the  other  aircraft. 

4)  Modified  LES  230-Y  Lateral/Dlrectlonal  Case 
a)  Direction  1  (Figure  C15) 

This  response  is  predominantly  roll  rate  over  all  times.  Beyond  4 
seconds,  roll  rate  and  roll  angle  components  are  roughly  equal,  while 
there  Is  slight  opposite  yaw  rate  beginning  around  2  seconds.  The 
sample  trajectories  (Figures  C16  through  Cl  9)  show  that  less  and  K>.- 
control  is  required  as  t^.  Increases.  Also,  note  that  the  only 
command  available  is  pure  aileron,  as  the  right  and  left  tips' 
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t i a jcc J ur ies  ,jii‘  mirror  images  of  each  other.  This  is  also  consistent 
with  H.FaS'IAB  symmetry  restrictions  as  well  as  the  stated  definition 
of  aileron  deflection. 

b)  Direction  2  (Figure  C20) 

This  response  direction  is  mostly  rolL  angle  except  at  1.5  seconds, 
when  opposite  rolL  rate  briefly  dominates,  then  quickly  diminishes 
until  J  seconds.  At  that  time,  opposite  roll  rate  again  grows  until 
it  is  the  secondary  component,  with  roll  angle  still  dominating. 

c)  Direction  3  (Figure  C21) 

This  direction  is  dominated  at  all  times  by  negative  yaw  rate.  At 
about  1  second,  a  minor  positive  sideslip  angle  component  begins.  At 
about  3  seconds,  another  minor  negative  roll  rate  component  begins. 
This  plot  shows  the  first  evidence  of  some  numerical  problems  in  the 
software  at  around  9  seconds,  as  all  vector  components  abruptly 
shift  values.  This  is  due  to  numerical  ill-conditioning  of  the  A  and 
B  matrix  combination  for  the  aircraft,  which  results  in  the  final  time 
Grammian  and  its  inverse  being  incorrectly  computed  by  the  existing 
library  matrix  subroutines.  A  specially-written  Grammian  inverse  sub¬ 
routine,  which  made  use  of  the  singular  value  decomposition  of  the 
matrix,  eased  the  numerical  problems  to  some  extent  but  did  not  remove 
them  completely. 

d)  Direction  U  (Figure  C22) 

This  most  difficult  response  is  always  dominated  by  sideslip,  with  a 
minor  yaw  rate  component  beginning  around  2  seconds.  Numerical 
problems  again  appear  around  9  seconds,  as  discussed  above.  The  very 
poor  controllability  of  this  direction  is  shown  in  Figures  C23  and 
C24  (2.0  seconds)  where  extremely  large  state  and  control  magnitudes 
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ari'  ri'<)u  i  i  ill .  Also,  Figures  C2b  <uul  F2/  show  numerical  problems, 
since  tlic  state  trajectories  never  reach  x^  exactly.  This  is  also 
attributed  to  the  particular  routines  used  to  compute  the  Graminian. 
e)  Cost  Plot  (Figure  C29) 

Cost  n ends  still  lollow  those  of  all  other  cases.  However,  the 
numerical  problem  al tecta  the  highest  two  cost  lines,  around  9 
seconds.  Clearly  this  aircraf t/i light  control  system  combination  is 
very  badly  behaved  in  the  lateral /direc L ional  mode, 
traits  Discussion 

The  response  characteristics  of  these  two  aircraf t/control  system 
combinations  all  exhibited  certain  trends: 

1)  Each  ranked  response  direction  appears  to  be  divided  into  three 
distinct  areas:  a)  the  very  short  final  time  range  (0  to  2  seconds) 
where  the  state  directions  remain  somewhat  constant;  b)  a  "transition 
region",  usually  from  final  times  of  2  to  5  seconds,  where  the  state 
variable  components  are  changing  from  their  short-time  to  long-time 
constant  values;  and  c)  the  long  final  time  region,  later  than  5 
seconds,  where  the  directions  are  again  fairly  constant  in  the  state 
space,  but  different  from  the  previous  directions.  These  three 
"bands"  may  be  physically  interpreted  as  the  dominance  of  the  short 
period/Dutch  roll  modes  and  phugoid/spiral  modes  in  the  first  and 
third  tj.  bands,  respectively,  and  the  transition  between  the  two  modes 
during  the  second  band. 

2)  It  is  very  difficult  for  the  flight  control  system  to  drive 
the  aircraft  along  these  "worst"  directions,  since  they  usually 
combine  conflicting  state  variable  motions,  e.g.,  positive  angle  of 
attack  along  with  negative  pitch  angle  and  rate,  as  shown  In  Figure  C9. 
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)  roiii  -i  film's  v  i  >'wpo  i  nl  ,  these  directions  usually  represent  highly 
unorthodox  .iiri  r.il  i  motions  tiiat  are  neither  required  nor  naturally 
ueour  during  normal  aircraft  use.  Figures  Cll  and  C13  again  show  that 
larger  tinal  times  require  less  control  magnitude,  probably  due 
to  the  aircraft's  dynamic  response  approaching  with  relatively 
little  extra  help  from  the  control  system. 

Direct  Aircraft  Comparison 

These  aireraf t /control  system  combinations  were  also  studied 
using  a  direct,  "head  to  head"  comparison.  Each  aircraft  is  "flown" 
to  the  same  arbitrary  set  of  state  vectors  for  both  the  longitudinal 
and  lateral  cases,  and  the  state  and  control  trajectories  are  compared 
to  see  which  aircraft  required  the  least  control  and  state  variable 
displacements.  Four  different  final  state  vectors  were  used  (see 
Table  V),  each  at  three  different  final  times  of  2.5,  5.0,  and  7.5 
seconds.  Certain  final  time  examples  show  results  typical  of  results 


Table  V 

Direct  Comparison  Final  State  Vectore 


VECTOR 

cf 

U (FT/SEC) 

BB 

^(DEG) 

h  (FT) 

LONG. 

7.5 

-5.0 

3.0 

0.25 

5.0 

25.0 

LONG,  1)2 

5.0 

-25.0 

8.0 

12.0 

8.0 

75.0 

VECTOR 

cf 

(3  (DEG) 

(DEG) 

-p  (DEG/SEC) 

V  (DEG/SEC) 

LAT.  til 

7.5 

0.0 

60.0 

25.0 

0.0 

LAT.  1)2 

2.5 

10.0 

5.0 

2.0 

5.0 

- - - 
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j  lu  .ill  three  times,  wiLhiti  a  given  final  state  vector,  and  these 
samples  art.  included  tor  all  tour  st.iLe  vectors  in  Appendix  I).  These 
typical  examples  are  now  discussed  in  detail. 

1)  Long  Ltndinal  Vector  1  (Figures  D1  through  D4) 

Overall,  the  LES  baseline  aircraft  is  by  far  the  most  controllable 
and  hence  requires  Lhe  least  control  energy  to  maneuver.  The  experi¬ 
mental  tip  controls  require  much  more  deflection  than  the  baseline  LES 
contro 1 s . 

2)  Longitudinal  Vector  2  (Figures  D5  through  D8) 

Again,  the  LES  baseline  aircraft  is  by  far  the  most  controllable  for 
the  least  control  movement.  Now  the  experimental  LES  aircraft  shows 
control  magnitudes  just  over  linearity  limits. 

3)  Lateral/Directional  Vector  1  (Figures  D9  through  D13) 

The  LES  baseline  aircraft  is  still  the  best,  with  the  experimental 
wlngtip  system  comparing  very  poorly  in  terms  of  control  deflection 
magnitudes.  This  is  caused  by  the  experimental  system's  having  only 
aileron  control,  with  any  state  space  movement  corresponding  to 
ordinary  yawing/rudder  movements  available  only  from  the  yaw  coupling 
response  due  to  excessive  aileron  control.  This  result  is  certainly 
consistent  with  expectations,  given  no  aerodynamic  yawing  control 
independent  of  the  roll  controls. 

4)  Lateral /Pi rectlonal  Vector  2  (Figures  D13  through  D16) 

The  LES  baseline  aircraft  remains  the  best  and  the  aileron-only 
experimental  aircraft  again  fares  very  poorly,  so  much  so  that  It 
violates  linearity  restrictions  for  control  magnitudes. 

Experimental  Control  System  Evaluation 

This  study  shows  that  the  experimental  wingtip  control  system 
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provides  only  nut  g  inn  1  longitudinal  amtrol  and  totally  inadequate 
later.il/cliroitiiMi.il  nmtrol.  This  problem  is  initially  predictable 
1 rom  the  small  magnitudes  of  the  VLEXSTAB-c.omputed  elevator  and 


aileron  control  power  derivatives,  C  and  C.  ,  as  shown  in  Table  A. Ill, 

% 


compared  with  the  corresponding  baseline  LES  control  powers.  The 
problem  is  also  quantitatively  demonstrated  in  the  direct  comparison 
tests  just  mentioned  above.  Due  to  the  FLEXSTAB  symmetry  restrictions, 
wingtip  rudder  power  has  not  been  accurately  modeled,  but  in  light  of 
the  wingtips'  correctly-modeled  yet  very  poor  elevator  and  aileron 
power,  the  rudder  power  Is  certainly  expected  to  be  equally  poor  with 
respect  to  the  baseline  LES  aircraft.  With  no  known  software  available 
to  model  nonsymmetrical  rudder  deflections,  there  is  no  method  to 
quantitatively  verify  this  expectation.  The  only  way  to  Improve  the 
yaw  control  capability  of  the  experimental  aircraft  so  that  FLEXSTAB 
could  correctly  model  its  control  characteristics  is  to  add  a  conven¬ 
tional  fin  and  rudder  on  the  aircraft's  centerline.  However,  this  is 
essentially  a  conventional  control  system  except  for  unusual  placement 
of  the  elevons.  Wind  tunnel  testing  appears  to  be  the  only  way  of 
gathering  data  for  any  asymmetrical  aircraft/control  deflection 
geometry. 

Method  Evaluation 

The  software  used  in  this  study  does  an  excellent  job  of  modeling 
linear  aerodynamics  of  an  arbitrary  (symmetric)  aircraft  geometry, 
then  translating  the  aerodynamics  into  a  state  space  linear  dynamic 
system  suitable  for  analysis  by  a  user-selected  linear  systems 
analysis  program  -  GRAMOP  In  the  case  of  this  study. 
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IIowl-vi- r ,  tin.  io  a  10  some  resit  irt  ions  which  must  be  observed.  The 
first,  FI.MXSTAK’s  symmetry  i  est  r  tc  L  ion  ,  has  been  described  already  so 
will  not  be  lurther  repeated  liere. 

Second,  the  linearity  restriction  of  all  software  must  be  under¬ 
stood.  Reference  1  and  Chapters  Ill  and  IV  describe  the  linearity 
assumptions  used  to  derive  the  three  programs  used.  Focusing  on  the 
GRAMOT  program,  it  is  necessary  to  understand  the  implications  of 
linearity  to  this  particular  optimal  control  problem  in  order  to 
correctly  interpret  its  solution.  The  constraint  equation  (34),  used 
in  the  optimal  control  analysis,  contains  matrices  A  and  B  which  are 
derived  using  these  1 inear ity /small  perturbation  assumptions: 

1)  linear  potential  theory  in  FLEXSTAB  aerodynamic  paneling; 

2)  small  perturbations  around  reference  flight  condition  to 
simplify  the  equations  of  motion.  Thus  the  matrices  A  and  B  are  a 
valid  physical  model  of  the  aircraft  only  for  those  aircraft  and 
control  motions  x(t)  and  u(t)  which  remain  within  the  region  of  small 
perturbations  with  respect  to  the  selected  equilibrium  flight  condi¬ 
tion.  Therefore,  the  only  physically  valid  GRAMOP-computed  optimal 
state  and  control  trajectories  are  those  which  also  remain  in  the 
region  of  small  perturbations  with  respect  to  initial  state  and 
control  vectors. 

Although  linearity  is  assumed  during  the  formulation  of  every 
software  algorithm  in  this  project,  several  cases  have  been  noted 
where  control  and/or  state  variable  magnitudes  have  violated  these 
perturbation  limits.  These  results  are  caused  by  the  combination  of 
two  effects  -  constant  A  and  B  matrix  elements  and  poor  controllability. 

First,  the  elements  of  both  the  A  and  B  matrices  are  constant 


41 


42 


uppropi  lately  limited  ran  ;>■  ol  o(  .  Figure  5C  shows  that  the- 
amstjnt-vd  1  no  uH'lliciein  assumption  erroneously  allows  the  algorithm 
to  extrapolate  the  stability  and  control  characteristics  beyond  their 
region  of  physical  validity. 

Second,  tire  CRAMOP  algorithm  makes  extensive  use  of  »  the 

diagonal  matrix  of  eigenvalues  of  W^(t).  Since  the  eigenvalues  have 
no  direct,  physical  state  space  meaning  in  terms  of  perturbation 
magnitudes,  they  are  only  limited  in  that  they  must  be  non-negative 
real  numbers.  But  this  matrix  and  its  inverse  appear  in  both  u*(t) 
and  aLso  x(t),  equations  (46)  and  (35).  As  described  before,  the 
relative  magnitudes  of  the  are  a  measure  of  controllability 

quality  along  each  of  the  orthogonal  eigenvector  directions  in  state 
space.  For  even  an  assumed  linear  system,  certain  ^ ^  may  be  so 

small  (i.e.,  those  eigenvector  directions  'f  may  be  so  poorly 

controllable)  that,  when  Inverted  in  equations  (46)  and  (35),  they 
drive  the  u(t)  and  x(t)  values  beyond  linearity  limits. 

This  combination  of  poor  controllability  and  extrapolated 
equations  of  motion  results  in  the  invalid,  nonlinear  output  data 
shown.  The  correct  interpretation  of  this  situation  is  that  GRAMOP 
cannot  always  compute  the  physically  correct  u*(t)  and  x(t)  outputs 
to  reach  every  possible  x(t^).  This  problem  exists  even  for  x(t^)’s 
which  are  themselves  within  small  perturbation  limits,  e.g.,  if  this 
final  state  lies  in  a  nearly  uncontrollable  subspace.  The  data  in 
such  cases  is  still  usable  qualitatively,  rather  than  quantitatively. 
Violation  of  small  perturbations  In  the  output  means  that  the  physical 
system  cannot  reach  the  given  x(t^)  when  constrained  within  only  the 
"linear  region  of  motion."  This  is  a  valid,  although  only  qualitative, 
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1L  is  certainly  possible  that 


measure  of  aircraft  ioiit  rol  lab  i  l  i  t  y . 
tltis  x(t  .)  could  be  reached  using  larger  perturbations  in  u*(t)  and 
x(t),  but  this  type  ot  motion  is  beyond  the  scope  of  both  FLKXSTAB  and 
GRAMOP .  In  order  to  acquire  valid  numerical  data  from  this  software, 
the  magnitude  of  the  "linearly  accessible"  x(t^.)  should  be  scaled  as 
small  as  necessary  until  the  largest  values  in  the  resultant  u*(t)  and 
x* (t)  trajectories  also  satisfy  the  linearity  constraints.  That 
scaled  x(tf)  then  represents  the  maximum-magnitude  of  the  "linearly 
controllable"  response  in  that  particular  subspace  of  (Rn) . 

The  third  restriction  is  that  all  solutions  to  this  optimal 
control  problem  depend  directly  on  the  user-input  control  weighting 
matrix,  R.  Even  If  all  constraint  matrices  A  and  B  and  all  other 
optimal  conditions  xf  and  t£  are  identical,  resulting  x*(t)  and 
u*(t)  trajectories  and  J*(u(t>)  will  differ  case-by-case  with  each 
different  R  matrix  used. 

Recall  that  the  R  matrix  form  adopted  in  this  study  is  diagonal 
with  non-negative  diagonal  elements.  This  choice  of  R  yields  a 
performance  functional  as  the  weighted  sum  of  squares 

J.suv,  *'**)** 

The  numerical  values  for  each  R  ^  directly  measure  the  relative 
control  energy  "penalty"  that  the  user  wishes  to  attach  to  each  u^. 
Each  Rti  establishes  a  deflection  magnitude  u^|  Is  re¬ 

garded  as  costing  one  "unit  of  equivalent  control  energy".  The 
quantity  u^'  -  u^/  4^  then  expresses  u^  as  a  "normalized"  u^',  in 
terms  of  "equivalent  control  energy"  units,  instead  of  physical  units 
(e.g.,  degrees  of  control  surface  deflection).  The  cost  functional 
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then  becomes 


(60) 


2 

Thus,  Rti  »  ,  or  the  square  of  the  arbitrary,  user-defined 

physical-to-equivalent  energy  units  conversion  factor.  The  larger 

becomes,  the  heavier  the  user  is  weighting  (penalizing)  the  use  of 
control  u^,  as  measured  by  J(u(t)).  R  allows  the  user  to  define 
equivalent  control  energies  for  controls  which  differ  physically. 

The  R  matrix  used  in  this  study  equated  the  control  energies  between 
angular  control  surface  deflections,  measured  in  degrees,  and 
throttle  movements,  measured  in  percent  of  aircraft  weight.  The 
author  arbitrarily  chose: 

1)  5  degree  control  surface  deflection  to  be  one  unit  of  control 
surface  energy 

2)  10X  of  aircraft  weight  to  be  one  unit  of  throttle  energy. 
Thus,  since  s  ( . 


r  o.o  t 

o.  o  1 

R  = 

L  0.0 

loo  J 

(61) 

rj,  1 

U1 

(62) 

For  the  lateral/directional  cases,  where  only  control  surfaces  (and 
no  throttle)  are  used,  R  was  left  as  the  Identity  matrix,  since  the 
same  ^  was  assumed  for  all  u^'s  in  these  cases. 
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VI.  Gone]  unions  and  Recommendations 

lurk  lu.siuns 

The  specific  conclusions  reached  during  this  study  are: 

1)  \  quantitative  measure  of  linear  system  controllability 
characteristics  has  been  developed  using  the  concept  of  control 
energy,  as  defined  by 

Jtcutn  =  ^  ^  sm  RwtTi  <ir 

0 

This  quadratic  functional  is  used  as  the  performance  index  in  the 

linear-constraint  optimization  problem;  find  the  u(t)  which  minimizes 

t 

=  i  I  u(t)Ru(riir 

0 

subject  to  the  constraint,  x  “  Ax  +  Bu,  with  x(0)  “  0  and  both 

t^  and  x(t^)  specified.  The  controllability  Grammian 

r  t  T 

=  j  #^)< lr 

0 

appears  in  both  u*(t)  and  x(t)  solutions  for  the  optimization  problem. 
When  W^ft)  is  expressed  by  its  unique  singular  value  decomposition 
form  within  the  u*(t)  and  x(t)  expressions,  it  gives  a  ranking  of  the 
controllability  quality,  i.e.,  control  energy  J(u(t)),  required  to 
drive  the  linear  system  along  each  of  the  n  unique,  time-varying, 
mutually  orthogonal  directions  in  state  space  described  in  Chapters 
I  and  IV. 

2)  Computer  software  program  GRAMOP  has  been  developed  which 
analyzes  any  input  linear  system  according  to  its  controllability 
quality  by  the  above  process  1).  Specific  program  options  include: 

a)  Calculate  and  display  the  time  histories  of  all  n 
orthogonal,  ranked  response  directions  at  any  given  final  time; 
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I>)  C.i  1 .  ii  l.i  l  i'  opilm.il  t  i  in  L  i  i>1  .mil  tor  respond  ing  state  variable 
response  u  jjci'li'i  ios  to  reach  any  input  linal  state  x(t^)  at  the 
spec i lied  i  ina  l  time. 

J)  CiRAMOi  has  been  saeeessl  ull  y  used  on  two  basic  aircraft/ 
flight  condition  combinations  which  were  compared  by  first  calculating 
longitudinal  and  lateral /direct ional  controllability  characteristics 
using  option  da  above.  Second,  optimal  control  and  state  response 
trajectories  were  then  calculated  to  achieve  selected  x(t^) 's  along 
each  ranked  response  direction,  using  option  2b  above.  See  Chapter  V 
for  details  of  these  results. 

4)  All  linearity  assumptions  inherent  in  the  constraint 
equation  (34)  must  be  carried  throughout  the  interpretation  of  any 
output  data.  In  these  aircraft  test  cases,  the  small  perturbation 
assumption  restricts  physically  valid  output  data  to  those  same  small 
perturbation  limits.  Any  "large"  output  data  is  correctly  inter¬ 
preted  as  showing  the  aircraft’s  inability  to  reach  the  desired 
x(tf)  when  using  only  "small"  control  inputs  and  state  variable 
responses.  Once  either  control  or  state  variable  output  magnitudes 
exceed  this  region,  the  linear  dynamics  model  is  an  invalid  repre¬ 
sentation  of  the  physical  aircraft  system.  In  such  cases,  the  mag¬ 
nitude  of  the  desired  x(t^)  should  be  scaled  down  until  linear/small 
perturbation  constraints  are  met  by  all  control  and  state  output 
data. 

5)  As  implemented  on  the  particular  LES  230-Y  aircraft,  the 
experimental  wlngtip  control  system  was  unsuccessful.  This  was 
proven  quantitatively  by: 

a)  very  low  FLEXSTAB-computed  values  for  control  power 
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tie i  iv.it  i  vfs ; 

Li)  poor  relative*  per l'ormauce  in  a  "head-to-head"  comparison 
with  the  baseline  aircratt; 

c)  fre«|uent  violation  ui  linearity  restrictions  throughout 
the  controllability  analysis.  The  control  system  symmetry  limitation  in 
FLKXSTAB  also  precluded  accurate  modeling  of  the  rudder  characteristics  of 
the  wine t i ps . 

Kecomme  ml  at  ions 

These  recommendations  center  around  changes  in  the  physical 
system  being  analyzed,  and  also  possible  modifications  to  the  CRAMOP 
algorithm. 

First,  this  study  looked  at  only  the  "raw"  stability  and  control 
derivatives  of  the  basic  aircraft  tested.  A  useful  extension  of  this 
analysis  would  be  to  modify  the  A,  B  matrix  representation  of  the 
aircraft  by  including  the  dynamics  of  additional  physical  control 
system  equipment,  such  as  Che  transfer  functions  for  servos,  stability 
augmentation  systems,  etc. 

Second,  different  cost  functionals  could  be  studied.  The 
present  J(u(t))  penalizes  only  the  magnitude  of  u(t),  and  allows  any 
u(t)  rates  with  no  penalty.  Better  aircraft-related  analysis  might 
result  from  including  some  control  rate  penalty.  Similarly,  penalties 
for  x(t)  and  x(t)  could  also  be  added.  This  would  make  the  cost 


Third,  other  software  changes  could  be  made  to  FLEXSTAB.  A 
routine  could  be  added  to  produce  linearized  state  equations.  A  most 


48 


.'Ii. 1 1 1  eng i ng  i  omputer/aerodynamic  project  would  he  to  replaire  Ei. EXSTAB '  s 


exist  ini;  nerodynam  i  c  paneling  code  with  other  more  advanced  paneling 
schemes . 
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Table  A. T 


Aircraft  Data  Specifications 


BOTH 

LF.S  -  230Y 
MODELS 

WING  AREA  (FT**2) 

670 

REF.  CHORD  (FT) 

25.55 

WING  SPAN  (FT) 

34.87 

AIRCRAFT  WEIGHT  (LB) 

36970 

MOMENTS  OF  INERTIA  (SLUG-FT**2) 

IXX 

25940 

IYY 

136800 

IZZ 

155500 

IXZ 

11700 
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Table  A. IT 


Trim  Parameters 


BOEING 

LES 

EXP'L 

LES 

MACH  NUMBER  (ND) 

0.80 

0.80 

ALTITUDE  (FT) 

35000 

35000 

VELOCITY  (FT/SEC) 

779 

779 

DYNAMIC  PRESSURE  <LB/FT**2) 

224 

224 

ROLL  ANGLE  (DEG) 

• 

0 

ROLL  RATE  (DEG/SEC) 

0 

0 

YAW  RATE  (DEG/SEC) 

0 

0 

NORMAL  LOAD  FACTOR  (ND) 

1.0 

1.0 

ANGLE  OF  ATTACK  (DEG) 

1.85 

2.11 

FLIGHT  PATH  ANGLE  (DEG) 

0 

0 

PITCH  ATTITUDE  (DEG) 

1.85 

2.11 

THRUST  (LB) 

5627 

2450 

ELEVATOR  ANGLE  (DEG) 

4.3 

35.7 

SIDESLIP  ANGLE  (DEG) 

0 

0 

AILERON  ANGLE  (DEG) 

0 

0 

RUDDER  ANGLE  (DEG) 

0 

0 
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Table  A. Ill 


KLEXSTAB  Stability  and  Control  Derivatives 


BOEING  LES 

EXP'L  LES 

cL  (ND) 

0.0906 

0.0661 

CDo  (ND) 

0.0231 

0.0046 

CMq  (ND) 

0.0320 

0.0386 

CLl  (ND) 

0.2464 

0.2464 

CDl  (ND) 

0.0375 

0.01633 

CMl  (ND) 

0.0 

0.0 

C,  (1/RAD) 

0.5446 

0.5722 

CDu  (l/RAD) 

0.0458 

0.0400 

Cj^  (l/RAD) 

-0.0272 

-0.0276 

CL^  (1/DEG) 

0.0377 

0.0332 

CD^  (1/DEG) 

0.0070 

0.0068 

CM<  (1/DEG) 

-0.0009 

0.0011 

CL  (l/RAD) 

2.1825 

1.8951 

CD  (l/RAD) 

0.1741 

0.1636 

CMq  <1/RAD> 

-0.7779 

-0.6509 

Cyp  (1/DEG) 

-0.0104 

-0.0025 

C#  (1/DEG) 

-0.0010 

0.0007 

Cn  (1/DEG) 

l> 

Cyp  (l/RAD) 

0.0030 

-0.0006 

-0.1572 

0.0931 

Cjp  (l/RAD) 

-0.1821 

-0.1276 

Cnp  (l/RAD) 

0.0860 

-0.0277 

Cyr  (l/RAD) 

0.7309 

0.0200 
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Table  A. TIT  (continued) 


_ _ _] 

BOEING  LES 

EXP  * L  LES 

C|  (l/UAD) 

0.1384 

1 

0.0488 

Cll(.  (  1  /RAD) 

-0.3679 

-0.0458 

C|.  .  (  L/RAD) 

<* 

0.7875 

1.1347 

l'u.  (1/HAD) 
oC 

0.0995 

0.1304 

CM.  ( 1 /RAD) 

<* 

-0.2344 

-0.3244 

Cy •  (I /RAO) 

V 

Ct.  (1/RAD) 

-0.2187 

-0.2205 

0.0033 

0.0018 

Cn ^  (L/RAD) 

CL^  (l /DEG) 

-0.0037 

-0.0104 

0.0196 

0.0031 

CD^  (1/DEG) 

0.0015 

0.0001 

Cmi  (1/DEG) 

E 

-0.0071 

-0.0011 

Cy  (1/DEG) 

0.0018 

-0.0016 

Cj  (1/DEG) 

cnc  (1/DEG) 

** 

0.0034 

0.0009 

-0.0011 

0.0005 

Cy^  (1/DEG) 

0.0044 

0.0002  * 

C i  (1/DEG) 

“R 

0.0009 

-0.0000  * 

Cn(  (1/DEG) 

R 

-0.0024 

-0.0001  * 

*  These  derivatives  correspond  to  strake  deflection 
and  are  ignored  in  any  subsequent  calculations. 


Append  is  11:  Sample  Output  Data  for  I.KS-2  30 Y 

This  is  a  collection  of  sample  output  data  generated  for  the 
Boeing  I.KS  2* 'JO—  V  using  Llie  original  flight  control  system  at  the 
flight  conditions  outlined  in  Chapter  V. 
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Figure  Bl  Baseline  LES  Longitudinal  'Ti(tf) 


Figure  B2  Baseline  IES  Longitudinal  Response  to  Reach  ^  Seconds) 
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Figure  B3  Baseline  LES  Longitudinal  Controls  to  Reach  ?t(tf  =  IS  Seconds) 
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Figure  B7  Baseline  IES  longitudinal  <?r(tf) 
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Figure  B8  Baseline  LES  Longitudinal  Response  to  Reach  *c  (tf  =  Seconds) 
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Figure  B9  Baseline  LES  Longitudinal  Controls  to  Reach  ^  (tf  =  IS  Seconds) 
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Figure  bio  Baseline  LES  Longitudinal  Response  to  Reach  ft  (tr  =  4.0  Seconds) 
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Figure  B13  Baseline  IES  L ate ra 1/Di  recti onal  ( t^) 
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Figure  bi4  Baseline  LES  Lateral/Oirectional  Response  to  Reach  ^x(tf  =  2.*r  Seconds ) 


Baseline  LES  lateral /Directional  Controls  to  Reach  v*  (tf  -2.2S  Seconds) 
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Figure  B17  8ase1ine  LES  Lateral/Directional  (tf) 
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Figure  B18  Baseline  LES  lateral/Directional  (tf) 
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Figure  bi9  Baseline  LES  Lateral/Directlonal  Response  to  Reach  v?4(tf  *  2.o  Seconds) 


SSKfT  LATERAL  CASE  CONTROL  INPUTS  VERSUS  TlftE 


Figure  b20  Baseline  LES  Lateral /Directional  Controls  to  Reach  ^  (tf  =  2.0  Seconds) 
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Figure  B21  Baseline  LES  Lateral/Di rectional  Response  to  Reach  ***  (tf  *  9°  Seconds) 
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Appendix  C:  Sample  Output  tor  Modified  l.ES  230-Y 

This  is  a  collection  of  sample  output  data  generated  for  the 
Boeing  l,KS  230-Y  using  the  proposed  new  flight  control  system  at  the 
flight  conditions  outlined  in  Chapter  V. 
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LEt  or  n  .8  36  RFT  UM01TUD1ML  CASE  EASIEST  REBrONSE  VECTOR  VS  TinE 


Figure  cl  Experimental  LES  longitudinal  (tf) 
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Figure  C2  Experimental  LES  Longitudinal  Response  to  Reach  vrx  (tf  =  f-F  Seconds) 
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Figure  C3  Experimental  LES  Longitudinal  Controls  to  Reach  (t^  =  Seconds) 
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Figure  c4  Experimental  LES  longitudinal  Response  to  Reach  (t*  =  4.o  Seconds) 
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Figure  C6  Experimental  LES  Longitudinal  trx  (tf) 
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Figure  C7  Experimental  LES  Longitudinal  v>  (tf) 
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Figure  C8  Experimental  LES  Longitudinal  (tf) 
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Figure  c9  Experimental  LES  Longitudinal  v?  (tf) 
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Figure  cio  Experimental  LES  Longitudinal  Response  to  Reach  u*  (tr  =  **°  Seconds) 
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Figure  ci3  Experimental  IES  Longitudinal  Controls  to  Reach  (tf  =  to  Seconds) 
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Figure  CIA  Log*n  of  Experimental  LES  Weighted  Control  Costs  for  All  Longitudinal  Basis  Vectors 


Figure  ci 5  Experimental  LES  lateral/Directional  (tf) 
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Figure  C16  Experimental  LES  Lateral/Directional  Response  to  Reach  O',  (tf  =  *•*  Seconds) 


Figure  ci7  Experimental  LES  Lateral /Directional  Controls  to  Reach  (tf  *  z,o  Seconds) 
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Figure  ci8  Experimental  LES  Lateral/Directional  Response  to  Reach  7,  (tf  *  4.0  Seconds) 
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Figure  C20  Experimental  Ltr>  Lateral /Directional  vv  (tj) 
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Figure  C21  Experimental  LES  Lateral/Directional  yj(tf) 
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Figure  c22  Experimental  LE!>  Lateral /Directional  (tj-) 
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Figure  C23  Experimental  LES  l.ateral/Directional  Response  to  Reach  (tf  =  z.o  Seconds) 
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Figure  C24  Experimental  LES  Lateral/Directional  Controls  to  Reach  <?4  (tf  *  z.o  Seconds) 


Figure  C26  Experimental  LES  lateral/Di rectlonal  Controls  to  Reach  ^  (tf  =  f.o  Seconds) 


figure  C27  Experimental  LES  Lateral/Directional  Response  to  Reach  O4  (t^  =  ^  Seconds) 


Figure  C28  Experimental  IES  lateral/Directional  Controls  to  Reach  (tf  *  ir  Seconds) 
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Figure  C29  Log)0  of  Experimental  LES  Weighted  Control  Costs  for  All  Lateral/Directional  Basis  Vectors 
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Append i x  1):  Samp le  Data  for  Direct  Aircraft  Comparison 

This  collect  ion  of  data  shows  a  direct  comparison  between  all 
three  aircraft.  They  were  each  flown  to  the  same  arbitrary  set  of 
two  longitudinal  and  two  lateral/directional  final  state  vectors,  aa 
outlined  in  Chapter  V. 
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Figure  di  Baseline  LES  Longitudinal  Response  to  Reach  Longitudinal  Vector! 


Figure  d2  Baseline  LES  Longitudinal  Controls  to  Reach  Longitudinal  Vector  I 
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Figure  d3  Experimental  IES  Longitudinal  Response  to  Reach  Longitudinal  Vector  1 
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Experimental  LES  Longitudinal  Controls  to  Reach  Longitudinal  Vector  1 
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Figure  D5  Baseline  LES  Longitudinal  Response  to  Reach  Longitudinal  Vector  l 
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Figure  r>6  Baseline  LES  longitudinal  Controls  to  Reach  Longitudinal  Vector  2 
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Figure  d7  Experimental  LES  Longitudinal  Response  to  Reach  Longitudinal  Vector  2. 


Figure  D8  Experimental  LES  Longitudinal  Controls  to  Reach  Longitudinal  Vector  2. 


UES  BMC  «  .B  SSKfT  LfiTCOW.  C09E  STfiTE  »E9fOKSE  ¥E»SUS  TlnE 


t®. 


B 


L®. 


•  1  1  1  "  1  ’ 
01*1 

"»  1  *■ 

0«' 

an 

r  .  T-.,T-r- 

DO*- 

iosoi  nawe  aneaoie 

X 

J  -L  -  (  . 

oo*m 

i  .  .  ,  J  . 

imt 

oi  •• 

•  1  *  * 

01*0- 

00  IS 

1030 1  3T0MU 

uuug 

OW1 

0 

1  1  1  i  ’  '  ' 
Oi*«t 

j  . . , . , 

0#  *12 

.1  . 

'  '  i  -  L  ’  ’  -  i 

MU- 

W«»- 

030/0301  SIMM 

no* 

V 

'  r  1 ) "  ’  1 

0|-| 

r  '  4  H  *  ' 
Of* 

■  •  '  t  •  ’  ■’!**" 

00*1 

*  |  M 

w»- 

00'«- 

1330/030)  3100  NBA 

X 

120 


Figure  09  8ase1ine  LES  Lateral/Directional  Response  to  Reach  Lateral  Vector  1 


Figure  dio  Baseline  LES  Lateral /Directional  Controls  to  Reach  Lateral  Vector  1 
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Figure  Dll  Experimental  LES  Lateral/DIrectional  Response  to  Reach  lateral  Vector  I 
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Figure  d12  Experimental  LES  Lateral/Directional  Controls  to  Reach  Lateral  Vector  1 
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Figure  di4  Baseline  LES  Lateral /Directional  Controls  to  Reach  Lateral  Vector  * 


Appendix  E: GRAMOP  Program  User's  Guide 


General 

This  guide  describes  the  detail  Instructions  for  using  the  GRAMOP 
computer  program  to  analyze  the  following  optimal  control  problem: 
find  the  optimal  control  history  u*(t)  which  drives  the  state  variable 
system  from  initial  state  x(0)  *  0  to  the  arbitrary  user-input  final 
state  x(t^),  and  which  also  minimizes  the  cost  functional 

n 

JcSfttt  =  z  1  aT<r)  Rstr)  Ir  (e-d 

o 

which  is  subject  to  the  state  space  constraint 

i  *  Ax  +  Bu  (E-2) 

Matrix  R  is  restricted  to  a  diagonal  matrix  with  non-negative  diagonal 
elements.  Details  of  the  theory  used  in  GRAMOP  may  be  found  in  Chapter 
IV  of  the  thesis,  "  Calculating  Aircraft  Controllability  Quality". 
GRAMOP  contains  two  basic  options: 

1)  TRAITS  option  -  calculates  the  controllability  characteristics 
or  quality  of  the  linear  system  throughout  the  n-dimensional  state 
space.  This  is  done  by  finding  the  n  time-varying,  mutually  orthogonal 
directions  in  (Rn)  made  up  of: 

1)  Vf^(t),  the  easiest  direction  of  all  to  control  in  (Rn); 

2)  VT  2^*)*  the  easiest  direction  of  all  to  control  in  (Rnl),  that 
particular  subspace  of  (Rn)  which  is  orthogonal  to  direction  1; 

1)  \r^(t),  the  easiest  to  control  direction  in  the  space 

(Rn  )  which  Is  orthogonal  to  direction  (i-1)  above; 
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n)  VT  (t),  the  most  difficult  direction  of  all  to  control  in  (Rn) 
n 

2)  XUHIST  option  -  calculates  both  the  state  history  x(t^,t^), 
and  the  optimal  control  history  u*(tg,tj)  which  drives  the  system  from 
Xq  =*  0  to  x(t^)  *  x^  along  the  state  trajectory  xft^.tj)  for  any  user- 
input  system,  xf  ,  and  t^.  This  option  can  be  used  to  investigate  some 
state  and  time  of  interest  to  the  user,  or  to  analyze  results  from  the 
TRAITS  option;  i.e.,  obtain  sample  control  and  response  histories 
which  move  the  system  along  a  specific  "ranked"  direction,  calculated 
in  the  TRAITS  option,  at  a  specific  final  time. 

Input  data  formats  are  as  follows: 


CARD  NO.  VARIABLE  NAME 

1  ISTOP 

DESCRIPTION  -  Plag  to  start 
CARD  NO.  VARIABLE  NAME 

2  N.M.I 

DESCRIPTION  -  N  -  number  of 
variables,  and  I  “  computer 


FORMAT 


List-directed  (LD) 

(1)  or  stop  (0)  GRAMOP  execution. 

FORMAT 

LD 

state  variables,  M  «  number  of  control 
hardware  flag:  (1)  *  INTERCOM,  (2)  -  batch. 


CARD  NO.  VARIABLE  NAME  FORMAT 

3,4  ID(1) ,  ID(2)  2A10 

DESCRIPTION  -  First  (second)  card  has  20  characters  for  first  (second) 
line  of  plot  title  box  of  all  plots 


CARD  MO.  VARIABLE  NAME  FORMAT 

5  NOUT  LD 

DESCRIPTION  -  Number  of  variables  x^  and/or  u^  to  be  plotted  or  printed, 
(must  be  less  than  or  equal  to  3). 
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FORMAT 


CARD  NO.  VARIABLE  NAME  FORMAT 

6,7  YTITLE  ( I , IM)  2A10 

YTITLE(1,20+IMI 

DESCRIPTION  -  First  (second)  card  has  20  characters  for  y-axis  print 
header/plot  Label  for  IMth  state  variable  (IMfc  cost  history).  There 
are  N  (see  card  2)  pairs  of  cards  6,7  -  a  pair  for  each  x . . 


CARD  NO. 


VARIABLE  NAME 


FORMAT 


8  YTITLE ( 1 , 10+IM)  2A10 

DESCRIPTION  -  contains  20  characters  of  y-axis  plot  label/print  header 
for  Tm^  control  variable.  There  are  M  (see  card  2)  cards  8,  one  for 
each  u^. 


CARD  NO. 


VARIABLE  NAME 


FORMAT 


9  IX(J),J-l,NOUT  LD 

DESCRIPTION  -  These  are  the  addresses/subscripts  J  of  all  x,  which  are 
to  be  printed  on  plotted.  (All  on  one  card  9)  J 


CARD  NO. 


VARIABLE  NAME 


FORMAT 


10  IU(J),J-l,MPLOT  LD 

DESCRIPTION  -  These  are  the  addresses/subscripts  J  of  all  u,  which  are 
to  be  printed  or  plotted.  (All  on  one  card  10)  ^ 


CARD  NO. 


VARIABLE  NAME 


FORMAT 


11  L1.L2.YY  LD 

DESCRIPTION  -  Matrix  A(B,R)  input  cards;  each  A  (B  ,  R  )  is  input 
as  i,J,A,,.  Each  card  except  the  last  ends  with  a^comma,  and  the 
last  card  r&r  each  matrix  ends  with  -  "0/".  There  is  one  set  of  cards 
11  for  each  of  the  three  matrices:  A,B,  and  R,  in  that  order. 


CARD  NO. 


VARIABLE  NAME 


FORMAT 


12  KTRAJ  LD 

DESCRIPTION  -  Output  option  flag;  (1)  *  plot  only,  (2)  «  print  only,  and 
(3)  “  both  print  and  plot. 


CARD  NO. 


VARIABLE  NAME 


FORMAT 


l  . 

I 


DESCRIPTION  -  This  is  the  XUHIST  routine  option  flag;  (1)  means  run 
XUHIST,  and  (0)  means  do  not  run  XUHIST  option. 

If  IXU  -  0,  go  to  card  19;  if  IXU  -  1,  the  following  cards  are  required: 
CARD  NO.  14  (TO, TF, NTPTS  LD) 

DESCRIPTION  -  TO  ■  user  specified  initial  time;  TF  *  user-specified 
final  time;  NTPTS  ■  number  of  evenly-spaced  "sample"  points  at  which 
output  data  is  calculated  in  the  interval  (T0,TF).  (Must  be  less  than 
100) 

CARD  NO.  15  x(TF)  LD 

DESCRIPTION  -  User-specified  final  state  vector  at  final  time  TF  (see 
card  14).  Same  format  as  card  11. 

CARD  NO.  16.17  (ID(3),ID<4)  LD 

DESCRIPTION  -  First  (second)  card  contains  20  characters  of  third  (fourth) 
lines  of  x(t)  plot  title  box.  Input  this  card  only  if  KTRAJ  (see  card  12) 
is  +  2. 

CARD  NO.  18  (NSCALE)  LD 

DESCRIPTION  -  The  x(t)  plot  y-axes  scaling  flag;  (0)  »  scale  all 
variables  Identically;  i.e.,  plot  all  state  variables  against  a  common  y- 
axis;  (1)  *  scale  all  variables  separately;  i.e,  each  variable  is  plotted 
against  its  own  scaled  y-axis. 

A  second  set  of  cards  16  -  18  now  follows,  in  identical  format,  except 
1D(3),ID(4)  are  now  the  titles  for  the  G*(t)  plot  title  box  lines  3  and 
4,  and  NSCALE  “  scaling  option  for  the  control  variables  history  plot. 

For  each  time  interval  (T0,TF)  and  final  state  (x(TF))  combination  desired, 
add  an  additional  set  of  input  cards  as  follows: 

a)  one  card  each  14  and  IS,  plus 

b)  two  sets  of  cards  16-18. 

To  exit  the  XUHIST  option,  input  a  final  card  16  as 
TO  »  0.  ,  TF  -  0.,  NTPTS  -  0 

CARD  NO.  19  (ITRAIT  LD) 

DESCRIPTION  -  This  is  the  TR/ITS  option  flag;  (1)  -  run  the  TRAITS 
option,  (0)  -  do  not  run  the  TRAITS  option. 

CARD  NO.  20  (TF1,TF2,NTF  LD) 

DESCRIPTION  -  These  Inputs  correspond  to  the  same  variables  as  those  on 
card  14. 

CARD  NO.  21  ( IDBUG  LD) 


ni 


■y 


DESCRIPTION  -  print  flag;  (1)  •  write  debug  print  onto  local 

file  name  TAPE7;  (0)  »  do  not  write  any  debug  print. 

Note:  Debug  print  includes  eigenvector  time  histories,  and  so  is  very 
desirable  to  have  for  later  XUHIST  use.  Remember  to  include  the  cards 
REWIND,  RERUN,  COPYSBF,  RERUN,  OUTPUT  in  the  JCL  statements  after  the 
GRAMOP  execution  card,  in  order  for  the  print  file  to  be  routed  to  the 
printer. 

CARD  NO.  22,  23  (ID(3),ID(A)  LD) 

DESCRIPTION  -  The  first  N(see  card  2)  pair  of  cards  contain  the  20 
characters  of  the  third  (fourth)  line  of  each  ranked  eigenvector  plot 
title  box;  the  last  pair  of  cards  contains  the  20  characters  for  the 
ranked  cost  plot  title  lines  3  and  A. 

A  new  set  of  cards  20  -  23  is  input  for  each  time  interval  (TF1.TF2) 
desired.  To  exit  TRAITS,  input  a  final  card  20  with  TF1  -  TF2  -  "0." 
and  NTF  -  "0". 

After  all  option  cards  for  the  TRAITS  option,  the  next  (or  last)  card 
is  another  card  1  with  ISTOP  -  0  to  end  the  GRAMOP  run,  or  with  ISTOP  * 
1  to  read  another  complete  set  of  GRAMOP  run  input  cards. 
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Sample  ORAMOr  Input  Heck 


ci$f?TITl  (..-•'»• 
CASE  T  ITL  '  (  f  n 

4 

ST  AT  £  VAH£r<  I  I 

oi  recti  ok  i  rorr 

STATE  VA>TABt'  * 

direction  2  <.or  r 

*f  AT c  VAFIAB1  r  ' 
DIRECT  I OK  1  COST 
STATE  VARIABLE  «* 
DIRECTION  4  COST 
CONTROL  1 
CONTROL  2 


\\ 

2 

•»»  •* 

i 

1 

-•9‘*5Elt  -f  t 

i 

? 

« 45 ?4e^-t  t 

i 

3 

,  33  7ieF~'l 

i 

4 

-•9<>f  441 ♦  C  j 

2 

1 

ft. 

2 

? 

0. 

2 

3 

♦  1  I  t:  tf  l 

2 

4 

.323E1P-M 

3 

1 

-.80437F+C1 

3 

2 

•  *>7  »2br;-f  3 

3 

3 

-.82'  491  ,G  1 

3 

4 

.r  3nr,t 

4 

1 

•  5i  474c  ^P l 

4 

2 

-.7  . 7  4t  E-t  4 

4 

3 

•  2<*B23E-«  ? 

4 

4 

-,23f  43C*€;. 

0/ 

1 

1 

•  17  64E-02 

1 

2 

•42454E-01 

1 

8* 

3 

c 

1 

•  4'578t*e2 

3 

2 

•6hZ7BE*C1 

•  284  7  3£*€ 1 
-•40450E*61 

1.  2  2  1.  !i> 


ri6:»IT.  39  2  1  •»  *174914  3  i 
STATE  VECTOR  PATH 
TO  XF  VS  TINE 

OPTIMAL  CONTROL 
INPUTS  VS  TIME 


i  •  > 


• » 


e.,o 

13«  «5b 


EASIEST  RESPONSE 
VECTOR  VS  TIME 
2N0  EASIEST  RESPONSE 
VECTOR  VS  TIME 
3RD  EASIEST  RESPONSE 
VECTOR  VS  TIME 
4T M  EASIEST  RESPONSE 
VECTOR  VS  TIME 
AL0G18  OE  EACH  OASIS 
COST  INDEX  VS  TF 

|» »#•»  • 


%  1  ,6ft tit  1/ 
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